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SUMMARY 


A method is presented for determining the effect of time lag in an 
automatic stabilization system on the lateral oscillatory stability of 
an airplane. The method is based on an analytical -graph! cal procedure. 
The critical time lag of the airplane -autopilot system is readily 
determined from the frequency-response analysis. 

The method is applied to a typical present-day airplane equipped 
with an automatic pilot sensitive to yawing acceleration and geared to 
the rudder so that rudder control is applied in proportion to the 
y® w ^^-S acceleration. The results calculated for this airplane -autopilot 
system by this method are compared with the airplane motions calculated 
by a step-by-step procedure. 


INTRODUCTION 


Recent calculations and flight tests of several airplanes designed 
for operation in the transonic speed range have indicated unsatisfactory 
damping of the lateral oscillation. The results presented in refer- 
ence 1 show that the oscillatory stability can be improved by the use 
of an automatic pilot. The calculations of reference 1, however, were 
made on the assumption of an idealized control system without lag. 
Reference 2 points out that lag of the type in which the amount of 
control applied at a given instant is assumed to be proportional to a 
deviation which existed at a fixed time previous to the given instant 


e‘ T s"b, 


can be represented mathematically by use of the lag operator 
where T 3 is nondimensional time lag based on the span and D^ is the 
differential operator. Lag of this type is generally referred to as 
time lag. Reference 2 suggests that for purposes of calculation the 
lag operator may be approximated by three terms of the Taylor's series 

Tu ** 


for 


s^b 


This approximation was used in some recent calculations 
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(reference 3) and the results were found to be erroneous and misleading. 
The purpose of this paper is to present a satisfactory method for 
determining the effect of time lag on the lateral oscillatory stability 
based on the exact expression of the lag operator rather than any 
approximation. Some recent analyses on the same problem, unknown to 
the authors at the time this problem was being analyzed, are presented 
in references 4 and 5.' 


SYMBOLS AND COEFFICIENTS 


0 

♦ 

P 

r 

V 

p 

V 

V 
P 

<1 

b 

S 

W 

m 

g 

^b 


angle of roll, radians 

single of yaw, radians 

angle of sideslip, radians (v/V) 

yawing angular velocity, radians per second (dty/dt) 

yawing angular acceleration, radians per second per second 
(d 2 t/dt 2 ) 

rolling angular velocity, radians per second (d0/dt) 
sideslip velocity along lateral axis, feet per second 
airspeed, feet per second 
mass density of air, slugs per cubic foot 

dynamic pressure, pounds per square foot ^ipV 2 'j 

wing span, feet 
wing area, square feet 
weight of airplane, pounds 
mass of airplane, slugs (W/g) 

acceleration due to gravity, feet per second per second 
relative -density factor (m/pSb) 
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rj inclination of principal longitudinal axis of airplane with 

respect to flight path, positive when principal axis is 
above flight path at nose, degrees 

7 angle of flight path to horizontal, axis, positive in a 

climb, degrees 


k Xo 

k Zo 

K *o 

K Zo 

K X 

K Z 


radius of gyration in roll about principal longitudinal axis, 
feet 

radius of gyration in yaw about principal vertical axis, feet 

nondimensional radius of gyration in roll about principal 
longitudinal axis ^kxQ^b j 

nondimensional radius of gyration in yaw about principal 
vertical axis ^kg 0 / /b^ 

nondimensional radius of gyration in roll about longitudinal 
stability axis ^|/k^ o ^cos^T) + K^^sin^T}) 

nondimensional radius of gyration in yaw about vertical 
stability axis (\Z k Zo 2cos2t 1 + K X 0 2 sin 2 r l) 


K. 


XZ 


C L 

c l 

^n 


nondimensional product -of -inertia parameter 


(K 


2 2 \ 

K Xq Jsin T] cos T} 


) 

trim lift coefficient j 

rolling-moment coefficient ( Rolli ^ b moment ) 
yawing -moment coefficient ^ ^ aw i n g^mome nt j 

/ Lateral forced 


lateral-force coefficient 


qS 


SC i 


Cz P = W 
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BCy 

=5T 


SC 
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■ , n T 


■'Hr 


A rb 

d 2V 


SC 

S 


n 

pb 

2V 


C 2 -*i 

P S** 
2V 


(3Cy 
' °2V 


c Y r = 


Scy 


^rb 

°2V 


SC 


Z 


' l r -srb 
3 2V 

Sc„ 


Cn s r = se; 


t 

s b 

Db 


1/2 


time, seconds 

nondimensional time parameter based on span (Vt/b) 
differential operator , 

period of oscillation, seconds 

time for amplitude of oscillation to damp to one-half its 
original value 
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deflection of control, radians 


CD 


real part of complex root of characteristic stability- 
equation 

angular frequency, radians per second 


CD- = — CD 

s v 


X = a + icD s 

T 


time lag between signal for control and its actual motion, 
seconds 


's-l' 


( T s) 


s'c 


Ka 


critical time lag 


maximum amplitude of acceleration in yaw produced by control 

- ) 
airplane/ 


deflection of unit amplitude 


< K a) s ■ ^ % 


amplitude of control-surface oscillation produced by auto- 
pilot in response to oscillation of airplane acceleration 

'l&r 


_r 


autopilot 


k s = ^ k 


K c =b 


( K c) s - r- 


phase angle, radians 
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* *• 

phase angle of lag of 6 r behind t when oscillating 
control surface forces airplane to oscillate, radians 

©q phase angle obtained from frequency response of autopilot 


ANALYSIS 


The investigation of the effect of time lag on the lateral 
oscillation may be conveniently divided into two parts: (a) determi- 

nation of the smallest time lag which would result in a neutrally 
stable oscillation, referred to as critical time lag, and (b) the 
effect of a given time lag on the lateral oscillatory stability. It 
is important to know some of the results of the analysis obtained in 
part (a) in order to facilitate the analysis presented in part (b). 

Part (a) is based on the frequency-response method of analysis (refer- 
ences 6 and 7) . This method affords a relatively simple means of 
determining the critical time lag of an automatic stabilization system 
and thereby of establishing the range of time lags for which the air- 
plane motion is stable. Part (b) treats the solution of a transcen- 
dental equation by means of an analytical-graphical procedure. The 
analysis is presented for sin airplane equipped with an automatic pilot 
sensitive to yawing acceleration and geared to the rudder so that 
rudder contfol is applied in proportion to the yawing acceleration, as 
suggested in reference 3* A similar analysis is applicable, however, 
to any automatic stabilization system with time lag. 

The equations of motion used are expressed in terms of the 
nondimens ional time parameter based on the span of the air- 
plane s^ = but the results of the calculations obtained in terms 

of sb have been converted from the nondimens ional time s^ to 

time t in seconds. Thus the discussion of the results and the figures 

included in the paper are given in terms of t in seconds. 


Determination of Critical Time Lag 

The critical time lag of a system is defined as the time lag that 
results in a neutrally stable or steady-state oscillation. The motion 
of the control 5 r and the airplane acceleration V when a critical 
time lag exists are shown in figure 1. This figure indicates that the 
relationship between the time lag t and the phase angle 0 between 
the motion of the control and the airplane acceleration can be expressed 
d 

as r = — . Since r can be expressed as a function of 0, the 
frequency-response method can be utilized to determine the critical 
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time lag for a combined airplane-autopilot system. This method 
requires that the frequency-response curves be calculated for the 
airplane and autopilot separately and the results analyzed to ascertain 
the conditions for neutral oscillatory stability of the airplane- 
autopilot system. 


Frequency-response curves for airplane .- The frequency-response 
curves for the airplane are obtained from the calculation of the steady- 
state motion of the airplane in response to a sinusoidal forcing 
function of unit amplitude (see reference 8). Thus if & r = sin o) s sb, 
the acceleration of the airplane is Db 2 ^ = (K a ) s sin(o> s Sb + 0 A ). The 
values of (K A ) s , known as the amplitude ratio, and 0 A are obtained 
over the desired range of angular frequencies by substituting ico s 
for Db in the expression for Db 2 ^/6r which is derived from the 
lateral equations of motion. The nondimensional lateral equations of 
motion, referred to stability axes, for a given control deflection are: 


( 2 ^K x 2 Db 2 - | C lp Db)0 + ^dbKxzV - | C Z r ^b) ^ ' C l^ = 0 
(^b^Z^b 2 - | 0^)0 + ( 2 ^K Z V - | C^Dbjt - C npP = C n5r & r 




KD 


("IcYp 1 ^ - Cl| 0 + ^2pbl>b " | Cy r Db " C L tan r)i|r + ( 2 HbI>b - c Yp)p = 0 

BC i BCv 

The derivatives = ^r- and Cy^ = are usually very small and 

therefore have been neglected in equations (l). 
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Hence, 



2 pbKx 2 Hb 2 - | C Zp Db 

0 . 


I>b 2 

2 hbKxZPb 2 - | Cj^Db 

c n 6 

u r 

-% 


- | Cy p Db - Cl 

0 

2 pbDb - Cy p 


5 r 

2p b Kx 2 Db 2 - \ C Zp Db 

2 ^b K XZ D b 2 - | C ^Db 



2 ^b K XZ I> b 2 - \ c n p D b 

2pbK Z 2 I>b 2 - | Cn r Db 



- \ CY p Db - Cl 

(2^b - | C Yr )l>b - Cl tan 7 

2|ibDb - Cy p 


( 2 ) 


After the numerator and denominator are expanded by the method of 

determinants, the expression for D^-,% / & r results in the ratio of two 
polynomials in D5. The substitution of ia) s for D5 in equation ( 2 ) 

. . id A 

gives a complex number A + iB which may be expressed as 


The amplitude ratio (%) s , which is equal to 
from the relation (Ka) s = \Ja‘ 
mined from the relation 6 


k a 


' 2 + B 2 . 


V' 


2’ 


can be determined 


The phase angle 0 A can be deter - 


. -IB 
y A = tan 


Frequency response of autopilot .- The frequency response of the 
autopilot is obtained from the equation for the control motion with 
time lag taken into account 

5 r = k s Db 2 \Ksb - t s ) (3a) 

where the term Db^Ksb - t s ) signifies the fact that the amount of 
control applied at a given instant is proportional to the acceleration 
at a fixed time t s previous to the given instant. This time-lag 

-t D 

effect can be expressed by the so-called lag operator e s D . Thus 
equation (3a) becomes 

S r = k s Db 2 ie" TsDb 


(3b) 
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Solving equation (3b) for 


Db^/Sr 


gives 


where 


Et! = Jl e TsDb = (Kc) £ e TsI>b 
5 r k s 


w 




6r 


autopilot 


Substituting ico s for Db in equation (4) results in an amplitude 
ratio (Kc) s = — , which is independent of frequency, and a phase 

"'kg . 

angle 9 C = TgO) s . 


Conditions for stability .- The necessary and sufficient conditions 
for any one of the oscillatory modes to be neutrally stable are that, 
at a particular frequency, the phase angles and amplitude ratios of the 
airplane and autopilot must be equal; that is, = 9q and (Ka) s = (Kc) s 


Even though these conditions are satisfied, the resultant motion of the 
airplane which is composed of all the individual modes of motion may be 
unstable since, as is pointed out subsequently, an additional unstable 
oscillatory mode may be present. In general, the condition which must 
be satisfied in order that all the oscillatory modes be stable is that. 


at each angular frequency where 0 a = ®C> the ratio 


(k a ) 


(Kc), 


— must be le 


ss 


than 1 at that frequency. The mathematical proof of this statement is 
given in reference 9* 


Illustrative example .- The foregoing method is applied to a typical 
present-day high-speed airplane having the characteristics presented in 
table I. The value of the control -gearing ratio k is arbitrarily 
assumed to be 0.0427. (This value of k = 0.0427 corresponds to a 
rudder deflection of 1° for a yawing acceleration of 23.4 D /sec/sec. ) 

The amplitude -ratio and phase-angle curves for the airplane, plotted as 

a function of angular frequency cu, where u> = c% Y., are presented as 

b 

solid-line curves in figures 2(a) and 2(b), respectively. As the 
frequency increases to infinity, Ka approaches a value of 15*98 
and (?a approaches it. The dashed line in figure 2(a), which is 
independent of frequency, is the amplitude ratio of the autopilot Kq. 

The phase-angle curves of the autopilot are straight lines with slopes 

equal to t, where T = - T s , and are shown as dashed lines in figure 2(b) 
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for several values of t. Since the phase angle remains between the 
range of 0 to 2it, 9q = tu> continues to repeat itself whenever tu 2 2n. 
To take account of this fact, is plotted as a series of parallel 
lines for each value of t. Figure 2(a) indicates that Ka = Kc 
at os = 3*8 and 03 = 8.5. The corresponding values of t 
where = at u) = 3-8 and 03 = 8.5. are t = I .63 and t = O. 38 , 
respectively. One of the oscillatory modes of motion is thus neutrally 
stable when t = I .63 and t = O. 38 . However, as mentioned previously, 
the motion of the airplane is neutrally stable only if all other 
oscillatory modes present are stable. This condition is satisfied 

for t = O. 38 , since for each value of 03 where 9& = 9r, the 
Ka ^ 

ratio — — <1. When T = I. 63 , one of the oscillatory modes is 

Kg 

neutrally stable but the system is unstable, because at cd=6, 


k a 

but ^>1. 


An analysis indicates that for values of Kc < 15*98, 


which is the limiting value of Ka, the system will be unstable for any 
infinitesimal time lag. The reason for the instability is that, for 
any infinitesimal time lag, the value of ©a is equal to 9q at some 

Ka 

high frequency where it can be shown that — > 1 since at the very 

Kc 

high frequencies 

KA _ (Ka^-^co + ak a 

K C (Ka)^- AKc 


where AKa and AKc are small incremental values. 


In order to verify the results predicted by the preceding analysis, 
motions of the airplane were calculated subsequent to an initial 
di sturbance of 5° in sideslip. The calculations involved a step-by- 
step procedure based on the Kutta 3/8 method (reference 10). The 
results of these calculations for t = O .38 and T = I .63 are pre- 
sented in figures 3 and 4, respectively. The solid-line curve in 
figure 3 was obtained by using a time increment of 0.095 and the motion 
is seen to be slightly unstable, whereas neutral stability is predicted 
by frequency-response analysis for t = O. 38 . An additional calculation 
was made by using a time increment of 0.0475, represented by the dashed- 
line curve in figure 3 , and, although the motion was still slightly 
unstable, the trend indicated by reducing the time increment was such 
as to make the oscillation more nearly neutrally stable. The airplane 
motion for t = I .63 is presented in figure 4 and, as was predicted, 
the motion is unstable. A neutrally stable oscillation was also 
predicted for this value of time lag but it is apparent from figure 4 
that the unstable mode influences the airplane motion more than the 
neutrally stable mode. 
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Effect of Time Lag on Lateral Oscillatory Stability 

Derivation of equations .- The nondimens ional equations of motion, 
referred to the stability axes, which include the effect of an auto- 
pilot applying rudder control in proportion to the yawing acceleration 
at time s b - t s > are obtained by combining equation (3b) with equa- 
tion (l). When 0 o e Xs k is substituted for 0, i|r 0 e^ & b for \j>, • 

and 3 0 e X& b for 3 in the resultant equation written in determinant 
form, X must be a root of the characteristic stability equation 

AX^ + BX3 + CX 2 + DX + E + k s e _T s X (a'x^ + B.'x 3 + C ’x 2 + d'x) = 0 (5) 

where A, B, C,D,E, A , B , C , and D* are functions of the mass 
and aerodynamic parameters of the airplane. The expressions for A, B, 
C, D, and E are given in reference 1 and 


A' = -Vb 2 Kx 2 C n5r 
B = (2p b K x 2 C Yp + n b C 2p ^C n5r 

c ' * (-s * s ‘V , f>) c S- 

D = 


The damping and period of the lateral oscillation in seconds are 
given by the expressions 


t 1/2 . | 


p = 2n b 
cp 6 V 




(6) 


where a and * co s are the real and imaginary parts of a complex root 
of equation (5). 


Determination of roots of transcendental equation .- The character- 
istic stability equation of this system (equation (5)) is seen to be a 
transcendental equation because a constant time lag t s in the automatic 

stabilization system is represented by the so-called lag operator e -Ts ^. 
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It is apparent that the complex roots of such an equation cannot be 
determined by conventional methods. A method of obtaining the complex 
roots of this transcendental equation for a particular value of r s 
is therefore presented. 

If a + io> s is substituted for X in equation ( 5 ) and the real 
and imaginary quantities are separated, two equations in a and cd s 
result: 


e" T s a = Fi(a,a> s ) 


(7a) 


e 


- T s a 


F 2 (a,o> s ) 


(7b) 


A simultaneous solution of equations (7a) and (7b) gives the desired 
values of a and oo s which satisfy equation (5). The method used to 
solve equations (7a) and (7b) simultaneously is basically a graphical 
procedure. For a series of values of o> s , the right-hand sides of 
equations (7a) and (7b) are plotted against a as illustrated in 
figure 5* The solid lines and dashed lines correspond to the func- 
tions F^(a,a) s ) and F 2 (a,a) s ), respectively. The variation of e -Tsa 
with a is also plotted in figure 5* The exact values of a and co s 

for which Fi(a,u> s ) = F 2 (a,<n s ) = e -Tsa are determined from a cross plot 
of the results of figure 5 as shown in figure 6. In this figure, 
with a as the abscissa and cjo s as the ordinate, are plotted the 

_T ; 

values of a and (d s which correspond to the intersection of the e s 
curve with the functions Fi(a,o) s ) and F 2 (a,o> s ). The solid curve in 

figure 6 thus satisfies the equation e” Tsa = Fi(a,o> s ) and the dashed 

-t a . . 

curve satisfies the equation e = F 2 (a,oo s ). The values of a 

and <d s at the intersection of these two curves therefore determine a 
root of the characteristic stability equation (equation (5))- The 
period and damping of the lateral oscillation are determined from 
equations (6) by using this root. 


An alternate method for the simultaneous solution of equations (7a) 


and (7b) is to substitute 1 - t a + 

s 


T 2 a 2 T 3a3 

— — — - — g — for e -T s a in each 


of these equations. The calculations involved become considerably less 
laborious when this approximation is made and results obtained by using 
this alternate method have been found to be in excellent agreement with 
results obtained from the method described previously. Figure 7 shows 
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that for | T s&| = 1> very close agreement is obtained between the exact 
_ t a. 

value, of e s and its approximation by three and four terms of the 
a 

series for e s . In actual practice the product T g a will almost 
invariably be much less than 1. When this substitution is made, both 
equations (7a) and (7b) become seventh degree in a for a given value 
of cu s . Although these equations are of high degree, no serious 
problem is presented since the values of a desired must be real and, 
in general, small; that is,, only one or two of the roots of these high- 
degree equations are of interest. In order to calculate the complex 
roots of equation (5) for a particular value of t s , the following 
procedure should be used. For a sequence of values of cd g , compute a 
from equations (7a) and (7b). The results obtained from each equation 
may then be plotted in a figure similar to figure 6 and the complex 
root of the characteristic stability equation determined from the 
intersection of the two resulting curves. If the value of the approx- 
imate series for e Tsa is in good agreement with the exact value 

of e -Tsa for the a determined from the intersection of the two 
curves, then the complex root obtained is valid. A point of inter- 
section for a value of a which would not give satisfactory agreement 

_ t o 

between the approximate series and e s might exist, however. If 
such be the case, the accurate point of intersection is readily 
ascertained. The estimated values of a and cd s represent the point 
of intersection of F^(a,co s ) or F 2 (a,o) s ) with the series approxi-r 

— T 8, 

mating e s . The desired point is the point of intersection 

of F^(a,o) s ) or F 2 (a,co s ) with e“ T s a . The first step is to evaluate 
_ r a 

the term e s for several values of a in the vicinity of the 
estimated point of intersection. The expressions F^(a,o> s ) 

and F 2 (a,cn s ) are then evaluated for values of co s slightly less than 

and greater than the estimated value of co s for several values of a 
in the range of the estimated point of intersection. Thus, the 
corrected curves of equations (7a) and (7b) are obtained and at their 
point of intersection, the accurate values of a and cn s are 
determined. 

A method for constructing curves of constant period and damping as 
a function of t s and k s is presented in the appendix. 

Range of o> s and a to be used in the determination of the com- 
plex roots of the characteristic stability equation .- In general, the 
analysis required to determine the damping of the airplane -autopilot 
system would be carried out for those values of t s that result in a 
stable system, since the purpose of equipping an airplane with an 
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autopilot is to stabilize, or increase the damping of the airplane 
motions. The range of t s for which the system is stable is deter- 
mined from the frequency-response curves of the airplane and autopilot 
as indicated in the section entitled "Determination of Critical Time 
Lag." Thus, if the system is stable, only negative values of a need 
be investigated. Also, if two stable oscillatory modes of motion exist, 
the least stable one is of greatest interest; that is, the complex root 
with the smallest real part is the one of most importance. 

The estimated range of values for (o a that should be used in the 
analysis to determine the effect of some particular value of time lag, 
located between t s = 0 and the critical time lag ( T s ) c , on the 
damping of the oscillation is obtained from the known values of the 
frequency of the oscillation for the cases where t s = 0 and t s = (t s ) c . 
The. ima ginary part of the complex root of equation ( 5 ) , which becomes 
a quartic equation when T s = 0, gives the value of 03 s for the case 
of t s = 0. The value of co s for T s = ( T s)c is determined from the 
analysis presented in the section entitled "Determination of Critical 
Time Lag." As mentioned in reference 11, the frequency of the oscilla- 
tion decreases as time lag increases; thus for a value of 0 < T s < ( T s ) c , 
the est ima ted values of co s should include frequencies greater than the 
value of at t s = (t s ) c and less than the value of oo 3 at T s = 0. 

Airplane motions for various values of t . - For purposes of com- 
parison, the motion of the airplane in sideslip subsequent to an 
initial displacement in sideslip of 5° was calculated, with the use of 
a step-by-step procedure, for t = 0, 0.1, 0.2, 0.25> 0.287> and O.38. 

The results are presented in figures 8(a) to 8(f). These figures 
indicate that as t increases from 0 to 0.2 the period increases 
slightly, whereas the damping is markedly improved. The frequency of 
the oscillation is about 03 = 3*7* However, as t continues to 
increase, the presence of a high-frequency oscillation is noted in the 
motion since this oscillation becomes less damped and the low-frequency 
oscillation becomes more heavily damped. (See figs. 8(d) and 8(e).) 

For t = O. 38, the motion is neutrally stable at the high frequency 
of oj = 8.5 but the low-frequency oscillation does not appear in this 
motion since it is very well damped. In figures 9(a) and 9(8), the 
period and damping of the lateral oscillation for several values of t , 
calculated by the method discussed in this paper, are compared with the 
period and damping readily obtained from the motion calculations shown 
in figures 8(a) to 8(f). The very good agreement between the results 
presented in figures 9(a) and 9(h) would probably be improved if the 
increment selected for the step-by-step calculations were reduced. The 
trends indicated by these results, however, are applicable only to the 
particular airplane -autopilot system considered in this example and may 
not be generalized to any other arbitrary airplane -autopilot system. 
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CONCLUDING REMARKS 


A method is presented for determining the effect of time lag in an 
automatic stabilization system on the lateral oscillatory stability of 
an airplane. The method is applied to a typical present-day airplane 
equipped with an automatic pilot sensitive to yawing acceleration and 
geared to the rudder so that rudder control is applied in proportion 
to the yawing acceleration. The results calculated for an airplane - 
autopilot system by the method described <are in good agreement with the 
airplane motions calculated by a step-by-step procedure. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 

Langley Air Force Base, Va., October 28, 19^9 
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APPENDIX 

A METHOD FOR CONSTRUCTING CURVES OF CONSTANT PERIOD 
AND DAMPING AS A FUNCTION OF t s AND k s 


If an over-all picture is desired of the effect of time lag for a 
given value of k s or the effect of varying k s for a given time lag, 
construction of curves of constant period and damping with the use of 
the following method is recommended. 

The characteristic stability equation of the airplane -autopilot 
system (equation (5)) may be rewritten in the form: 


AX^ + BX.3 + CX 2 + DX + E i. _ -T s^' ( /\-\ \ 

-(a'jJ* ♦ B’x3 + C V ♦ D'l) = M (M) 

If X = a + ico s is substituted in each side of equation (Al), the 
condition that X be a root of the characteristic equation is that the 
complex number Aq + iBi obtained from the left-hand side must be equal 
to the one obtained from the right-hand side Ag + iB2. The 

quantities Aq + iBq and A 2 + iB2 may be represented by the 

expressions Rqe 101 and R2e 102 , respectively. Therefore, this require- 
ment is equivalent to saying that Rq = R 2 and 0q = @2 if ^ is to 
be a root of equation (Al). If X = a + ia> s is substituted in the 
right-hand side the following expressions result: 

i0 

k s e~ Tsa (cos - i sin T s a> s ) = Aq + i®2 = ^2 e ^ 


where 


R 2 = k s e- T s a 


and 


_q -sin t s (d s 

0 O = tan — — = -t-od = 2 n - t ca 

2 cos TgOJg s s s b 
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Therefore, if \ = a + iui s is substituted in the left -hand side 

i0i 

and A^ + iB^ = R^e x is obtained, the value of r required to 
make 0q = 02 can be determined. Since T s is therefore determined 
and a is fixed, the value of k s necessary to make Rq = R 2 can be 
calculated. Thus for these values of k s - and t s , \ = a + ico s is a 
root of the transcendental stability equation (equation (5 ) ) • 

For a given value of a the analysis may be made throughout the 
range of co s and the corresponding values of k s and t s determined. 

A curve representing this value of a may then be plotted as a function 
of k s and t s . This procedure is repeated for a sequence of values 
of a and the corresponding curves in the k s ,T s plane are plotted. 

Each point on a curve of constant a represents a particular value 
of o) s . Curves of constant frequency may therefore be plotted by 
drawing a curve through the given value of oj s on each one of the 
a curves. The values of a and o> s are converted to TWg and P 
by equations (6). The final result would consist of curves of constant 
damping T ^/2 and curves of constant period P in the k g , x s plane. 

Thus the effect of time lag on the lateral oscillations for any value 
of the- gearing ratio k s , or the effect of varying k s for any value 
of time lag, may be ascertained. For purposes of illustration several 
lines of constant Tq/2 and P in the k , t plane were calculated for 
the typical present-day airplane described in table I and are shown in 
figure 10. 

A curve can be plotted in the k s ,T s plane which divides the 
quadrant into a satisfactory and an unsatisfactory region according to 
any prescribed relationship between the period and damping of the lateral 
oscillation. In order to calculate this curve, several values of a 
and cd s that exactly satisfy the criterion should be selected and 
substituted in equation (Al). The combination of k s and t s is then 
obtained for each set of values of a and co s and the desired curve 
is plotted in the k s ,Ts plane. 
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TABLE I 

STABILITY DERIVATIVES AND MASS CHARACTERISTICS OF A 
TYPICAL PRESENT-DAY AIRPLANE 


W/S, lb/ft2 

S, £t 2 . . 

b, ft 

p, slugs/ft3 
V, ft/sec . 
7, deg . . 
Ct i • « • • 


*4 ■ ■ 

K Z 2 . . 
Kxz • • 

deg 

per 

CZ r , per 


per 
Cn r , per 


Cy p , per 
c Y r > P er 


radian . 
radian . 
radian . 
radian . 
radian . 
radian . 



per radian 


Cnp> per radian . 
C per radian . 
C n g > P er radian 


. . 65 
. . 130 
. . 28 
0.00089 
. . 797 
. . 0 
. 0.23 
. 80.7 
0.00967 
0.0513 
-0.001^5 
. - 2.0 
. - 0.40 

. 0.08 

- 0.0155 


- 0.40 
. . 0 

. . 0 


- 1.0 


, 0.25 
-0.126 


-0.163 
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(h) Phase angle. 


Figure 2.- Frequency-response curves of the typical 
present-day airplane and the assumed autopilot. 
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Figure 4.— Effect of t =1.63 on the motion of the airplane. 
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(c) T=0.20. 


(f) r* 0.38 (sameasf/g. 3). 


Figure 8.— Effect of time lag on the airplane motion in sideslip. 






28 


NACA TN 2005 



0 .10 .ao .30 40 

Tj sec 

(a) Damping. 


3.0 


2.0 


R sec 


1.0 


0 



0 JO .ao .30 .40 

Tj sec 


(h) Period. 

Figure 9* — Effect of time lag on the period and damping of the lateral 

oscillation. 










